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Abstract 



An idea which has been around in general relativity for more than forty years is 
that in the approach to a big bang singularity solutions of the Einstein equations 
can be approximated by the Kasner map, which describes a succession of Kasner 
epochs. This is already a highly non-trivial statement in the spatially homogeneous 
case. There the Einstein equations reduce to ordinary differential equations and it 
becomes a statement that the solutions of the Einstein equations can be approx- 
imated by heteroclinic chains of the corresponding dynamical system. For a long 
time progress on proving a statement of this kind rigorously was very slow but 
recently there has been new progress in this area, particularly in the case of the 
vacuum Einstein equations. In this paper we generalize some of these results to the 
Einstein-Maxwell equations. It turns out that this requires new techniques since 
certain eigenvalues are in a less favourable configuration in the case with a mag- 
netic field. The difficulties which arise in that case are overcome by using the fact 
that the dynamical system of interest is of geometrical origin and thus has useful 
invariant manifolds. 
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1 Introduction 

The fundamental equations of general relativity are the Einstein equations, possibly cou- 
pled to other equations describing the dynamics of the matter which generates the gravi- 
tational field. With a suitable choice of physical units the equations are 

Rot/3 — = T a p. (1.1) 

The unknowns in these equations are the spacetime metric g a p and the matter fields. 
R a p is the Ricci tensor of the Lorentzian metric g a p and R its trace. T a $ is the energy- 
momentum tensor. In this paper we are mainly concerned with the Einstein vacuum 
equations, where T a p = 0, and the Einstein-Maxwell equations. In the latter case the 
source of the gravitational field is an electromagnetic field F a p and the energy-momentum 
tensor is given by 

T aP = FJF^-^F^F^g^. (1.2) 

The electromagnetic field tensor is antisymmetric (F a p = —Fp a ) and satisfies the source- 
free Maxwell equations 

V a F aP = 0, V a F M + V 7 F Q/3 + V p F ia = 0. (1.3) 

It is well known that solutions of the Einstein equations generally develop singular- 
ities. In particular, there are solutions relevant to cosmology in which the singularity 
corresponds to the big bang. Belinskii, Khalatnikov and Lifshitz (hereafter abbreviated 
to BKL) developed a heuristic picture of the singularities in cosmological solutions of the 
Einstein equations. In this context they introduced a map of the circle to itself which we 
refer to as the Kasner map. It will be defined precisely below. They suggested that it 
provides a model for oscillations of the geometry in the approach to the singularity. For 
the original work see |BKL70] and |BKL82j . A modern discussion of these ideas can be 
found in [HU09J. An important idea in the BKL work is that spatially inhomogeneous 
solutions of the Einstein equations can be approximated by spatially homogeneous so- 
lutions near the singularity. Since from a mathematical point of view the dynamics of 
spatially homogeneous solutions is still far from understood it is natural at the present 
time to concentrate on understanding classes of spatially homogeneous solutions. This is 
the strategy we adopt in what follows. 
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A long-standing question in mathematical cosmology is to relate the Kasner map 
to the dynamics of actual solutions of the Einstein equations, possibly with matter. An 
important recent advance in this field is the paper [LHWGllJ where a relation of this kind 
was established in a special case. These results concern solutions of the vacuum Einstein 
equations of Bianchi types VIII and IX. They complement earlier results of Ringstrom 
|Rin00j . |Rin01j by providing a more detailed description of the dynamics of the approach 
to the singularity in certain cases. The work of Ringstrom on vacuum spacetimes was 
preceded by results of Weaver [WcaOOj on solutions of the Einstein-Maxwell equations of 
Bianchi type VIo using a dynamical system introduced in [LKW95] . The aim of this paper 
is to extend the results of |LHWGlT] to this case of the Einstein-Maxwell equations. There 
is other recent work on this question in the vacuum case |BeglO , jRTIO] but these papers 
use very different techniques from those which we will apply to the Einstein-Maxwell case 
and for this reason they will not be discussed further here. 

In the next section the necessary background and the fundamental equations needed 
in the paper are introduced. The most important similarities and differences between the 
models with magnetic fields considered in what follows and the vacuum models which had 
previously been analysed are explained. The third section contains the main theorem and 
an exposition of the strategy of its proof. The central result is the existence of unstable 
manifolds of codimension one for some heteroclinic chains. To prove this it is necessary 
to obtain estimates for a solution during its passages close to the Kasner circle and for its 
behaviour between passages. This is done in sections [4] and [5j respectively. A central idea 
of the paper and one which is a major step beyond what was achieved in the vacuum case 
is the use of a specially constructed Riemannian metric to measure the distance between 
the heteroclinic chains and the approximating smooth solutions. The last section discusses 
future extensions of this research and interesting open problems. 



2 The basic set-up 

Spatially homogeneous spacetimes are those solutions of the Einstein-matter equations 
where there is an action of a Lie group G by isometries of g a p with three-dimensional 
spacelike orbits which leaves the matter fields invariant. The cases where the isotropy 
group is discrete can be classified according to the Lie algebra of G. It is common in 
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general relativity to use the terminology due to Bianchi, who introduced types I to IX. It 
is also common to distinguish between two subsets of these types known as Class A and 
Class B. In what follows we will only be concerned with Class A models. More information 
on this subject can be found in [WE97] or |Ren08| . 

The analyses of vacuum spacetimes mentioned above are based on the well-known 
Wainwright-Hsu system |WH89] . This is a system of ordinary differential equations for 
five variables (£+, Ni,N 2 , N 3 ) which are subject to one constraint. It includes all the 
Bianchi models of Class A (i.e. types I, II, VIo, VIIo,VIII and IX). The system is defined on 
a smooth hypersurface in R 5 . An analogous system for Bianchi spacetimes of type VIo with 
a magnetic field was introduced in [LKW95J . It is also defined on a smooth hypersurface 
in R 5 and it includes solutions of types I and II with a magnetic field. The variables 
are called S + , £_, N + , N_, H. The first two variables can be identified with the variables 
of the same name in the vacuum case since they have the same geometrical meaning 
in both cases. The variables N + and N_ correspond in a similar way to certain linear 
combinations of N 2 and N 3 . More specifically, N+ = §(JV 2 + A^) and AL = &(N 2 - N 3 ). 
The variable H corresponds to the magnetic field. 

The dynamical system is 

E' + = -2A^(l + £+) + §# 2 (2-£ + ), 

£'_ = -(2N 2 + §# 2 )£_ -2N + N_, 

N' + = (2£+(l + £+) + 2£ 2 + §if 2 )AT + + 6£_AL, (2.1) 

N'_ = (2S + (l + S + ) + 2S 2 + §F 2 )AL + 2£_AL, 

H' = -(£+(2 - £+) - S 2 + N 2 )H. 

The prime denotes a derivative with respect to a time variable r which tends to — oo 
as the singularity is approached. These equations are taken from [WcaOOj. They arise 
as a special case of the equations for models with a magnetic field and a perfect fluid 
derived in [LKW95J by setting the fluid density to zero. (Here a magnetic field means an 
electromagnetic field satisfying the condition that F^rtr = 0, where n a is the unit normal 
vector to the group orbits.) Solutions are considered which satisfy the condition 

T? + + £ 2 + Nl + -H 2 = 1. (2.2) 

This condition follows from the Einstein equations and is preserved by the evolution 
equations for (£ + , A^., AL, H) just defined. The inequalities AL > 0, N 2 < 3N 2 and 
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H > are assumed. These are also preserved by the evolution and define the region which 
corresponds to Bianchi type VIo solutions with non-zero magnetic field. Setting H = 
while maintaining the other two inequalities gives a representation of the vacuum solutions 
of Bianchi type VI . Setting iV_ = or N + = V3N_ gives two different representations 
of solutions of type II with a magnetic field. Setting other combinations to zero leads to 
vacuum solutions of type II, solutions of type I with a magnetic field and vacuum solutions 
of type I (the Kasner solutions). Note the invariant subspaces {N 2 = 0} and {N 3 = 0} 
which also appear in the Bianchi system with perfect fluid. In fact the invariant subspaces 



{iV2 = 0}, {N$ = 0}, {H = 0} will play a crucial role in our analysis, see (3.4). 

The circle defined by + S 2 = 1 consists of stationary points. Each one of them 
corresponds to a Kasner solution and so this set is called the Kasner circle. There are three 
families of heteroclinic orbits between points on the Kasner circle whose projections to the 
(£+, S_)-plane are straight lines. Two of these families correspond to vacuum solutions 
of Bianchi type II and occur in both the vacuum case and the case with magnetic field. In 
the vacuum case there is a third family related to these two by symmetries of the system. 
In the case where a magnetic field is included the two sets of Bianchi type II vacuum 
solutions are complemented by a family of Bianchi type I solutions with magnetic field. 
The projections of the latter to the (£+, £_)-plane are identical to those of the third family 
of Bianchi type II solutions in the vacuum case. There is thus a natural correspondence 
between heteroclinic chains consisting of Bianchi type II solutions in the vacuum case and 
heteroclinic chains in the case with a magnetic field which include orbits corresponding 
to both solutions of the vacuum Einstein equations of Bianchi type II and solutions of the 
Einstein-Maxwell equations of Bianchi type I. In the vacuum case there is a heteroclinic 
cycle consisting of three orbits and it is the central example considered in |LHWGlT] . 
The projections of the orbits making up this cycle to the (£+, S_)-plane are related by 
rotations by multiples of \. By what has already been said, there is a corresponding 



heteroclinic cycle in the system of [LKW95J. See also figure 2.1 
The Kasner solutions can be written in the explicit form 



- dt 2 + t 2pi dx 2 + t 2p2 dy 2 + t 2p3 dz 2 . (2.3) 



With a suitable choice of ordering the Kasner exponents pi, p 2 , p 3 are related to the 
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Figure 2.1: Kasner circle and period 3 cycle; Bianchi type II vacuum families in black, 
Bianchi type I family with magnetic field in blue; Kasner intervals /Q bounded by Taub 
points Ti of tangencies and antipodal points Qi. 

variables £+ and £_ by 

Pl = |(1 -2E+), 

p 2 = |(1 + E + + V3E_), (2.4) 

p 3 = |(1 + E + -V^S_). 

The eigenvalues of the linearisation of the Wainwright-Hsu system at a Kasner solu- 
tion are (6pi,6p 2 ,6p 3 ) (|WH89), p. 1425). For the system describing solutions of the 
Einstein-Maxwell equations of Bianchi type VIo introduced in |LKW95] the eigenvalues 
are (3pi, 6p 2 , 6p 3 ) ( |LKW95j . (4.1)). The Taub points of the Kasner circle are defined by 
the condition that the Kasner exponents are (1, 0, 0) or a permutation thereof. All points 
of the Kasner circle other than the Taub points have a one-dimensional stable manifold. 
For pi < the stable manifold is defined by a solution with non- vanishing magnetic field 
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and the corresponding eigenvalue is 3pi < 0. For P2 < or p% < the stable manifold is 
defined by a vacuum solution. 

In order to have an overview of the relative sizes of the different eigenvalues it is 
useful to introduce the Kasner parameter u G [l,oo] which is defined implicitly by the 
relation (cf. |HU09j ) 

-u 2 (l+u) 2 . . 

VWiV* = 7z - - 2vr 2.5 

Then the Kasner exponents arranged in ascending order are given by 

P 1 ~ l+u+u 2 ' 

= irfe, (2-6) 

~ u(l+u) 
" 3 l+u+u 2 • 

Note that p2 > —pi with equality only when u = oo. On the other hand y < —p\ and 
Y > —p\. The Kasner map is defined by u t— > u — 1 for m > 2 and « 4 (u- l)^ 1 for 
u < 2. The central example in 1.1 1 Ij is a set of heteroclinic orbits which form a cycle 
of order three. In fact the value of u corresponding to that example is invariant under 
the Kasner map. It solves the equation u 2 — u — 1 = and is the golden ratio |(1 + \/5)- 
The values of the Kasner exponents at the vertex of this cycle where p\ <pi < Pz are 

Pl = |(i -V§), 

P2 = I (2-7) 
p 3 = 1(1 + a/5). 

For each value of the Kasner parameter u in the interval (1, oo) there are six points 
on the Kasner circle where u takes that value. Removing the Taub points and their 
antipodal points Qi from the Kasner circle leaves a union of six intervals Ki, 1 < i < 6. 
They will be numbered as follows. Let K\ be the region where p\ < and P2 > P3- Then 
number the others consecutively while moving anticlockwise along the Kasner circle, see 



figure 2.1 



In order to assess the stability of a heteroclinic cycle it is important to examine the 
eigenvalues of the linearisation of the system at the vertices. In both the vacuum case 
and the case with magnetic field there is one negative eigenvalue — p,\ and two positive 
eigenvalues /X2, /^3- Without loss of generality the labelling can be chosen so that /i2 < /i3- 
Then from what has been stated above it can be seen that in the vacuum case the 
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inequalities —pii < ^2 < ^3 hold at any point of the Kasner circle except Tj and Qi. 
Call this the first linearisation condition. That this is true is one of the most important 
hypotheses of the main theorem of |LHWGlT] . On the other hand this condition can 
fail in the case with magnetic field. It fails precisely when the Kasner exponent p\ is 
intermediate in size between p 2 and P3, i.e. when the point of the Kasner circle lies in 
the one of the sets K 2 and K 5 . For then the eigenvalue 3pi is smaller in magnitude 
than that of the negative eigenvalue, while the other positive eigenvalue is not. Then we 
have the situation that u. 2 < —fix < fis- The situation that the eigenvalue 3pi does not 
correspond to the eigendirection tangent to the heteroclinic orbit incoming towards the 
past is covered by the theorems in this paper. Call this the second linearisation condition. 
What is common to the first and second linearisation conditions is that the eigenvalue 
corresponding to the heteroclinic orbit incoming towards the past is larger in modulus that 
that corresponding to the heteroclinic orbit outgoing towards the past. In the example of 
the 3-cycle at least one of the two linearisation conditions just introduced holds at each 



of the the vertices. See figure 2.1, the second eigenvalue condition holds in the intervals 
K.2 and /C5, whereas the first eigenvalue condition holds in the remaining intervals. 

The linearisation conditions are not in themselves enough to make the theorems in 
this paper work. Additional geometrical information is required. This is the existence of 
a certain invariant manifold. It is tangent to the space spanned by the vectors tangent 
to the stable manifold and the centre manifold and the eigenvector corresponding to the 
largest eigenvalue. For a general dynamical system there is no reason why a manifold of 
this kind should exist. In the example of a Bianchi model of type VIo with magnetic field 
a manifold of this kind is defined by the vacuum solutions of type VI or the solutions of 
type II with magnetic field. 



3 Main result and sketch of proof 



We shall prove the following result on the dynamics of the Bianchi model of type VIq (2.1 



2.2) with magnetic field. 



Theorem 3.1 The period 3 heteroclinic cycle given by (2.1) possesses a local codimen- 
sion-one unstable manifold. In other words, system (2.1, 2.2) admits a codimension-one 
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manifold, locally close to the heteroclinic cycle, of initial conditions whose backward tra- 
jectories converge to the heteroclinic cycle. The manifold is locally Lipschitz continuous in 
the open complement of the boundaries N\ = 0, N2 = 0, H = 0. It is Lipschitz continuous 
in every closed cone intersecting these boundaries only in the heteroclinic cycle. 



The proof will only use certain properties of the particular structure of the Bianchi 



system (2.1) and can be sketched as follows. 



Step 1: local passage, section |4j In a neighbourhood of the equilibria of the 



heteroclinic cycle, i.e. close to the Kasner circle, the Bianchi system (2.1, 2.2) with 
reversed time direction can be smoothly transformed to a vector field 



x = f(x), x = (x u , x ss , x s , x c ) e R x E x R N x R 



that satisfies the following properties: 



Conditions 3.2 (loc-i) There is a straight line of equilibria, 



/(0,0,0,z c ) = 0. 



(3.i; 



(loc-ii) The heteroclinic orbits of the original system correspond to the x ss - and 
(loc-iii) The linearisation at the origin has the almost diagonal form 



( 



Df(0) 



flu 



\ 



"/•is 



-/isl 



* * 



-/^sN 

* J 



(3.2) 



with /tu,/t ss ,/t s i, • • • ,/i S N > 0. 

(loc-iv) The eigenvalue corresponding to the incoming direction is stronger than the 
eigenvalue corresponding to the outgoing direction, 



/t u //t SS < 1- 



(3.3) 
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(loc-v) The codimension-one subspaces 

{x u = 0}, {x ss = 0}, {x sl = 0}, . . . , {x sN = 0} (3.4) 
are invariant under the flow. 



Note that, for the model (2.1, 2.2), we have N — 1. However the generalization 
N > 1 is needed for applications discussed at the end of this section and in section [6j 

The last two properties are the most crucial for the proof. The invariance of the 
codimension-one subspaces is non-generic for systems admitting (loc-i)-(loc-iv) and is a 
very strong constraint on the system. 

Under the above assumptions the local passage map, \l/ loc : £ m — > S out , from an 
in-section E in = {x ss = e} to an out-section S out = {x n = e}, for sufficiently small 
s, is a Lipschitz-continuous map with arbitrarily little change of the x c -component and 
arbitrarily strong contraction transverse to the x c -component. Unfortunately, this only 
holds true with respect to a non-Euclidean metric and represents one of the main diffi- 
culties of our investigation. Indeed, for /i s i < /i u even a linear vector field / never gives 
rise to a Lipschitz-continuous map \l/ loc with respect to the Euclidean metric. The case 
/i u < /i ss ,/i s i, . . ./i S N, on the other hand, does yield a Lipschitz-continuous map ^ loc and 
has been treated in |LHWGlT] . 

Step 2: global excursion, section [5j Close to the heteroclinic chain, by smooth 
dependence on initial conditions, the trajectories follow the heteroclinic orbit from the 
out-section of a local passage to the in-section of the next local passage. This map, ^f° h ■ 
E£ ut — > S 1 fe n +1 , is a uniformly bounded diffeomorphism. In particular, any deformation 
imposed by i]>s lob in directions transverse to x c will turn out to be dominated by the 
strong contraction of the local passage map \l/ loc . In x c -direction, however, we gain an 
expansion given by the Kasner map. Thus, the global excursion ^f oh : S£ ut — > given 



by the Bianchi system (2.1, 2.2) with reversed time direction satisfies 



Conditions 3.3 (glob-i) ^f oh maps the origin of ££ ut to the origin of S 1 fe n +1 and (local 
neighbourhoods of of) the invariant subspaces {x ss = 0}, {x s i = 0}, . . . , {x s n = 0} 
onto {x u = 0}, {x s i = 0}, . . . , {x s n = 0} (in arbitrary order). 

(glob-ii) ^jf ob isaC 2 Diffeomorphism. The bounds ||Z^^| lob || , || J D(^ lob )- 1 ||, \\D 2 ^f oh \\, 
||D 2 ( 1 I'| lob ) _1 || < M are independent of k. 
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(glob-iii) It uniformly expands in x c -direction at the boundary. In other words, ^f oh : 
{x° s ut = 0, i™' = 0} — > {x™ = 0, = 0} is Lipschitz continuous and its inverse 
(^| lob |{£out =0iX out =0 }) _1 has Lipschitz constant less than L < 1, independent of k. 

Step 3: graph transform, section [5j Combining local passage and global excursion 
yields maps from each in-section to the next, \1/ = ^| lob ° \l/}? c : — > with uniform 

cone conditions. A standard graph-transform technique now yields the claimed invariant 
manifold as a fixed point in the space of Lipschitz-continuous graphs x k = x^(x^,x^) in 
S™. For completeness of presentation, we will give the necessary arguments in section [5] 



In fact, steps 1-3 prove a much more general theorem than 3.1 that is: 



Theorem 3.4 Let a C A vector field and a chain of heteroclinic orbits h k (t), 

lim h k -\t) =p k = lim h k (t), k e IN. 

t— >oo t— >— 00 

be given. Assume that locally near p k assumptions (loc-i)-(loc-v) hold and that along each 
h k assumptions (glob-i)- (glob-iii) hold, with constants a,L,M independent ofk. 

Then there exists a local codimension-one stable manifold to the heteroclinic chain, 
i.e. a codimension-one manifold of initial conditions following the heteroclinic chain and 
converging to it. 

The heteroclinic chain itself is contained in the boundary of this manifold. The man- 
ifold is locally Lipschitz continuous in the open complement of the invariant subspaces 
(loc-v). The manifold is uniformly Lipschitz continuous in every closed cone intersecting 
the invariant subspaces only in the heteroclinic chain itself. 

This theorem covers not only the period 3 cycle of the Bianchi VIo system with 
magnetic field. In fact, it applies to every heteroclinic chain in the Bianchi VIo system 
with magnetic field that does not accumulate at any Taub point (required for uniformity 
of bounds) and such that the chain does not contain heteroclinic orbits of the magnetic 



family to points in the domain K2 U K§. See also figure 2A_ It also applies to every 
heteroclinic chain of the Bianchi A (VIII and IX) system without magnetic field but 
with ideal fluid as investigated in [LHWGllJ, as long as it does not accumulate at any 
Taub point. The proof given here completes the arguments sketched in the discussion 
section of |LHWGiT] and relaxes the constraint on the matter model required there. All 
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matter models that yield positive eigenvalues of the linearisation at the Kasner circle in 
the non-vacuum direction are included, due to the relaxed eigenvalue condition. 



4 Local passage near a line of equilibria 



In this section we study the passage of trajectories under a general flow near a line of 



equilibria with eigenvalue constraint (3.3) and invariant subspaces (3.4) consistent with 



the Kasner circle in the Bianchi VIo system with magnetic field. We will collect estimates 
on expansion and contraction rates to establish Lipschitz properties of the local map 
between sections to a reference orbit given by the passage near the line of equilibria, see 



theorem 4.8 at the end of this section. Compared to |LH WGlT] (section 3), we assume the 
relaxed eigenvalue condition (3.3) without any constraint on . . . ,/i s n- This requires 



the use of a non-Euclidean metric (4.21, 4.22). 
Consider a C k vector field, k > 4, 



x 



X 



(x n ,x ss ,x s ,x c ) G R x R x r x R, 



(4.1) 



that satisfies conditions |3.2| in a neighbourhood of the origin. Due to the invariant sub- 
tle 

\ 



spaces (3.4), the form of the linearisation (3.2) holds locally all along the line of equilibria, 

/ 



D f (0,0,0, x c 



-/^l(x c ) 



V 



-/isN(£ c ; 
* 



o 



(4.2) 



The stable and unstable manifolds as well as the strong stable foliation of the stable 
manifold are C k and can be flattened, see e.g. |SSTC98] . Theorem 5.8. By a C k change 
of coordinates the stable / strong stable / unstable manifolds to the equilibria locally 



coincide with the respective eigenspaces. In particular, and in addition to (3.4), the 
following stable and unstable fibres become invariant: 



W u (x c ) = {x ss = 0,x s = 0, x c fixed}, 
W s (x c ) = {x u = 0,x c fixed}. 



(4.3) 
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out 




Figure 4.1: Local passage ^ loc : S in ->■ S out . 



Note that in the Bianchi system, W n (x c ) coincides with the outgoing heteroclinic orbit 
attached to the equilibrium (0, 0, 0, x c ). 



Due to (4.3), the linearisation becomes diagonal, 



Df {0,0,0, x c 



-/4i(x c ) 



-HsN\X c ) 

J 



(4.4) 



Our aim is to study a local map from an in-section E in = {x ss = e} to an out-section 
j-jout _ _ £ | £ Qr Xg;Xc ~ see fig ur e 4.1 This corresponds to the passage near the 



Kasner circle in the Bianchi system in backwards time direction. (We reversed the time 
direction to obtain a well defined local map.) 

We rescale the system to 



x = A(x c )x + eg(x), 



(4.5) 



with e arbitrarily fixed and g at least quadratic in (x u ,x ss ,x s ). Due to the invariant 
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subspaces (3.4) and ( |4.3 ), the vector field takes the form 



L...JV, 



(4.6) 



-fi s i(x c )x s e + eg s e(x)x si , 

e (g CBB {x)x Bg + Ylt=i9c S i{x)x s ^ x u , 

with C fc_1 -functions g u , g ss ,g S £, vanishing along the line of equilibria, and C fc_2 -functions 
g C ss,g CS £- In particular 

\gu(x)\,\g ss (x)\,\g sl (x)\,...,\g sN (x)\ < Cmax(|x u |, |x ss |, |z fll |, . . . , |x sN |), 
\g C ss(x)\,\g csl (x)\,...,\g csN (x)\ < C, 

for some constant C > independent of e and igW, where IA is some local neighbourhood 
of the origin. Similarly, all derivatives of g u , g ss , g S £, g css , g CS £ are bounded by C for x G U. 
We choose 

U = (-2,2) JV+3 . (4.8) 
All further estimates will use this rescaled system ( |4.5 ) with flattened invariant manifolds 



(4.3) in the local neighbourhood 14. They will be valid for all e < £q and suitably chosen 



Eq. In the original system (4.1), e bounds the size of the neighbourhood of the origin in 



which this local analysis is valid. 

Proposition 4.1 Let 

/i u :=/i u (0), -/i ss : = 



-/*ss(0), 



1 N, 



be the eigenvalues of (4-4) a ^ the origin. Then for all < a < 1 there exists an e > 



such that for all e < Eq in (4-5) and x G 14 

Hu{x c ) u. ss (x c ) /jL s i(x c ) 



a < 



/' 



fJ- S N[Xc) 



< a' 



(4.9) 



Proof. Due to the invariant subspaces (3.4, 4.3), the linearisation of the system at 



eigenvalues close to the origin remains diagonal, and the eigenvalues depend differentiably 



on x c , For the rescaled system (|4.5|) with small £0 this provides bounds in 14: Indeed, there 

< eC. (4.10) 



exists a constant C > independent of eo, e, such that 















1 


dXc ^ss( x c) 


1 


dt^( X c) 



[X] 
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The scalar function 9{x) := fi u {^u{x c ) + £9vl(%)) 1 is therefore C k 1 and close to 1. 
The vector field 



x ' = 6(x)f(x) 



:/(*) 



fx u (x c ) + eg u (x) 

has the same trajectories as the original vector field and all previous considerations remain 
valid. Thus we can assume, without loss of generality, that 6(x) = 1 in U, i.e. 



g u (x) = 0. 



(4.11) 



At this step we have made use of the fact that the origin possesses exactly one unstable 
eigenvalue. The vector field to consider then has the form 



X n /^u^-u 

x s i = -/j si (x c )x se + eg se (x)x s i, 
x c 



1,...N, 



(4.12) 



£ (g css (x)x ss + Y,eLi 9csi{x)x s ^j x n . 



Lemma 4.2 For all < a < 1 there exists an e > such that for all e < Eq, x(0) G U 
and t > 0, as long as x(t) remains in U under the flow to the vector field (4- Iffy , we can 
estimate: 



\x c (t) - x. 



xjt) = exp(/x u t) x u (0), 
x ss (t) e [exp (-i// ss t) ,exp (-0;^^*)] x ss (0), 
Xst(t) e [exp (-^JJ-stt) ,exp(-afiset)] x^O), 
2eC I 1 



< 



a 



/4s fJ'VL 



(4.13) 
(4.14) 

1,...,N, (4.15) 
(4.16) 



i/ere, C is the uniform (in x and e) bound from (4-7) 



Proof. The unstable component (4.13 ) is given directly by the vector field. The estimates 



of the stable components (4.14, 4.15) follow from the vector field and the uniform bounds 



(4.9, 4.7). Indeed for arbitrary < a < 1 and e < Eq small enough, we have 



■-/4s - eC, -afi ss + eC 
a 



Thus for arbitrary < a < 1 we find suitable a < a < 1 and Eq small enough such that 

1 



3 s 



a 



• ^ s 
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Integration yields the claim. Bounds on x S £ are obtained analogously. The centre com- 



ponent (4.16) is then estimated by plugging (4.13, 4.14, 4.15) into the vector field (4.12) 



and integrating: 
\x c (t) -x c (0)\ 



N 



9 CSS 

(s))x ss (s)x ll (s) + ^ 9cse(x{s))x si (s)x u (s) ds 
1=1 

< eC exp(/i u - q/i sa ) |x ss (0)3: u (0) | + ^ exp(-q/i s i) \x si (0)x u (s) \ ds 

Jo i=1 

nt N 

< 2sC / exp(^ u - 5>fi ss ) \x n (0) | + exp(-o> s£ ) |x s£ (0) | ds 

Jo i=1 

ft n 

< 2sC exp(-a(/i ss - /i u ))|x u (0)| + ^exp(-a/i S £)|x s ^(0)| ds 

Jo i=1 

The last inequality needs a slight adjustment of a ~^ a and uses the eigenvalue condition 



(3.3). Indeed, for all < a < 1, we find a suitable < a < 1 with < a(/x ss — /i u ) < 
aUss - A* u - cx 



The local map 



(V n „in i: 



^out out „OUt\ 



\l/ loc (a; in x in x inN ) 



(4.17) 



,,m 

■ y ss 



is given by the first intersection of the solution of (4.12) to the initial value (x^,x^ 
1, xf-, x™) with the out-section {x u = 1}. See figure 4.1 The local map \l/ loc is well-defined 
on the in-section 

S in = { («>r,4 n ) I *S = 1, < a£ < 1, HxMl < 1, |41 < 1 } (4-18) 



see lemma 4.3 below. The singular points in the intersection of the stable manifold of the 
equilibrium line with the in-section are mapped to the respective points in the intersection 
of the unstable manifold of the equilibrium line with the out-section: 



^ (-^u ~ J^c ) — (-^ss i X s i X c ) -~ (0)0, X c ). 



(4.19) 



Note that there is no drift in x c at the boundary due to the invariant fibres (4.3) 



Lemma 4.3 There exists an Eq > such that for all e < Eq and x(0) = x m in the in- 



section S in , see (4-18), the trajectory x(t) under the flow to the vector field (4-12) remains 
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inU as long as \x n \ < 1, i.e. all along the passage defining the local map \l/ loc . 
time t loc is given by 



t 



loc 



il,' 

/i u I X v 



The passage 
(4.20) 



Proof. We choose e smaller than 



i 



2(N+1)C 



amines - /i u , /i sl , . . . , /x bN }), see lemma 



4.2 



Then trajectories starting in E in cannot leave U unless x u becomes larger than 1, see 



(4.14, 4.15, 4.16). Furthermore, (4.13) ensures that x u must grow beyond 1. Thus every 
trajectory starting in E in intersects the out-section E out = {x u = 1} before leaving U. 



Setting x u (t loc ) = 1 in (4.13) determines the passage time t loc . 



Corollary 4.4 The local map \l/ loc (^.i7, 1^.19), i.e. the local passage on the closed in- 
sertion E* including the singular boundary {s™ = 0}, is continuous. For all < a < 1 
there exists an Eq > such that for all e < sq the following estimates hold 



|™out „in 

l x c x c 



< £ c (\x™\ + \x™\ + ■ ■ ■ + \x%\) , 



x, 



out I 

si I 



<- |™in|a/i ss //i u — 1 1 ™in I 

— l A u I l A u |) 

^ \^.in\afj, s i/ fj, u |™in| 

— I x u I l x s£h 



with C independent of e and x m . 

Thus the drift along the line of equilibria is arbitrarily small and the distance from 
the orbit to the union of the stable and unstable manifolds shrinks arbitrarily fast, close 
to the critical orbit. 



Proof. The estimates follow directly from lemma 4^ applied to the local passage time 
given by lemma 4.3 They also establish continuity of the local map ^ loc at the singular 
boundary {s™ = 0}. Note < x™ < 1 on the in-section. Note further that a/i ss /yU u — 1 > 
for a chosen close enough to 1. ixi 



Unfortunately, there is no hope to obtain Lipschitz estimates for the local map \l/ 



loc 



with respect to the standard metric. Even for g = 0, the linear vector field (4.5) yields a 
non-Lipschitz local passage, for fi s e < fi u . 

To obtain Lipschitz bounds for the local map \l> loc , we have to introduce a non- 



Euclidean metric on the in- and out-sections. We define the Riemannian metrics 

N 



dst 



1 3^1,8 1| 2 



da; 



= 1 



X 



(4.21) 



at. 
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on the in-section E m and 

II II 2 N II II 2 

dsl = dx l + dx s 2 £ + dxl (4.22) 

X ss e=1 X si 

on the out-section S out . We denoted the Euclidean norms ||£ u ,s||2 = II { x m x s) III = x u + 
J2t=i x lti \\ x ss,s\\i - \\( x ss, x s )\\l = x 2 s + Y^=i x lt Tne distance in £ m , S out is then given 
by the length of the shortest connecting paths and denoted by dist*. On fibres {x c fixed} 
we define the metric analogously. 

Let us discuss the new metric in the cone x UjS = (x u ,x s ) G [0, oo) N+1 in E m , ignoring 
the x c -direction that remains unchanged. The metric becomes singular along the invariant 
boundaries {x u = 0}, {x s e = 0}, see (3.4). Inside the open cone (0, oo) N+1 , the new metric 
ds* is locally equivalent to the Euclidean metric ds, 

ds 2 < dsl < (Ar + l) ^!^)' 1 ^ 1 '---' 1 ^ 1 ^ ^, (4.23) 

mm{|x u |, \x Bl \, . . . , |x sN |} 

and thus induces the same topology. The origin can be included. In fact the distance of 

any point to the origin is bounded by 

IKslh < dist*(0,x u , s ) < (iV + l) 3 / 2 ||x UiS || 2 . (4.24) 

(The upper bound can easily obtained by connecting the origin to x u ^ s with a piecewise 
linear path along the space diagonals with respect to suitable coordinate directions.) 
Every curve in the open cone hitting the boundary away from the origin has infinite 
length. 

In particular, the new metric is uniformly equivalent to the Euclidean metric in any 
closed cone that has finite, nonzero angle to the boundaries, i.e. 

{ (x n , B ,x c ) | ||xu, s ||2 > cmax{|x u |, |x s i|, . . . , |x sN |} }, c> 1. 

Thus Lipschitz estimates with respect to the new metric carry over to the Euclidean 
metric. 

We denote the in- and out-sections without the singular boundaries but with the 
origin by 

S m . = E m n((0,oo) w+1 U{0})xR) 

= { (xn,s, %c) e S in I x UjS = or x u x sl ■ ■ ■ x sN ^ }, 
S out . = s out n ((0, oo) N+1 U {0}) x E) 

= { (^ss,s5 x c ) G S out I x SSjS = or x ss x s i ■ ■ ■ x s n 7^ }. 
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Corollary 4.4 yields Lipschitz continuity of the local passage \& at the origin with 



respect to the Euclidean metric, and by (4.24) also with respect to the new metrics (4.21 



4.22). To obtain Lipschitz estimates away from the invariant boundaries with respect to 



the new metric, we consider the linearisation of the vector field (4.12) along a trajectory 
x(t) from the in- to the out-section to obtain bounds on D^/ loc . 

We start with a tangent vector 5 m = (8™, 8™ = 0, 8™, <5j. n ) of unit length with respect 



to the metric (4.21 ) at a point x 
1 = Il5 in ll? = 



x 



in 112 



\X 



U,sll2 I rin|2 



in 1 2 



G S" 

N 



x 



in 112 



X 



u,sll2irin|2 , |c 
ail 



in 1 2 
c I ' 



(4.26) 



First, we project 8 m along the vector field / into the hyperplane {8 n = 0}, as this remains 
invariant under the linearised flow and corresponds to the out-section. The projected 
vector 



8(0) 



8 in - 



8]] 



[i u X u 



(4.27) 



thus represents our initial condition to the linearised flow 

8 SS = (-u ss (x c ) + eg ss (x) + e(d Xas g ss (x))x S8 ) 8 SS 

+ e Y.e=i( d x Be 9ss(x))x ss 8 s e 

+ (-u' ss (x c ) + e(d Xc g SB (x))x ss ) 8 C , 
8 s e = (-^ s e(x c ) + eg a (x) + e(d Xsi g si (x))x ae ) 8 s e 

+ e(d Xas g s e(x))x s i8 BS 

+ £ 'Ei^t(9 XJ g Bt (x))x ae 6 B i 

+ (-/Xrf(a: c ) + e(d Xc g se (x))x se ) 8 C , 
8 C = e[d x c 9css(x)x ss x n 4" i d Xc g cs i(x)x S £X u j 8 C 

+ £ (g css (x)x u + d Xss g css (x)x ss x u + J2f=i d Xas g cs i(x)x se x v ) 8 S 

Here we already dropped the w-component. 



(4.28) 



U S £. 



Lemma 4.5 Let a unit tangent vector 8 m to x m = (x™ B ,x^) G S m with respect to the 



metric (4-21) be given. The projection 8(0) of 8™ along the vector field (4-12) into the 
plane {8 n = 0} is estimated by 

\8 SS (0)\ < c\kj 2 \ 

\8J0)\ < Clbl nl| - J| ' 



'u,sll2 



X 



s£l? 



\$c(0)-8?\ < eC, 
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with a constant C independent of x m , 5 m , and e < Eq, provided Eq is chosen small enough. 



Proof. Apply (4.27) to (4.26) and use the bounds (4.7) on the nonlinear terms of the 



vector field (4.12) 



Indeed, we find 



W = 5 m -j£kf(* m ) 



( s in 


\ 




5 in 


= 










We" 


/ 





(-/i ss (4 n )+^ ss (x in )) 1 
(-u, se (x?)+Eg s£ (x in ))x i » 



Immediately, we have 5™(0) = 0. For the other components we again use the uniform 



bounds (4.7) on the nonlinearity g, the bounds (4.9) on the eigenvalues, and the bounds 



(|4.26|) on the components of S m . We obtain for arbitrary < a < 1, if £ < £o is chosen 

1 



small enough: 



< 



IC 


1 




™in 


a 



/4s < 



"Mu IK n sll2' 



for the component transverse to the in-section, 



I*. 



< \%\ + 



K 


1 


/4i 


™in 


a 



fist\x^\ < 1 + 



the 



si I 



for each of the remaining N stable components, and 









™in 1 
x u 1 



1 



eC(N + 1)\x™\ < e—C 



< e—C, 



A'u IF u ,bI|2 

for the centre component. An obvious choice of a new constant C yields all claimed 



[X] 



estimates. 

Lemma 4.6 For all < a < 1 there exists an Eq > such that for all e < Eq, the 



linearised flow (4-28) can be estimated: 



|<U' < -a/i ss \8 ss \ +eC\x ss \ ( \8 C \ + J2t=i IM ; • 



\b~ s e\' < -afi s e\8se\ + £C\x si \ (\S c \ + \8 ss \ + J2fye \^d\ 



\8 C \ < eC (Jx ss | + Yli=i \ x si\) \ x n\%\ + eC\x u \ [ \5 SS \ + J2e=i \5 se \ , . 
Here C is a constant independent of x m , 5 m , and e < Eq. 
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Proof. Use the bounds (4.7) on the nonlinear terms of the vector field (4.12) and the 
bounds (4.10) on the derivatives of the eigenvalues. Note that igW = [— 2,2] N+3 . This 
immediately yields the claimed estimates. M 

Lemma 4.7 For all < a < 1 there exists an 6q > such that for all e < Eq the 
following statement holds: Let a trajectory x(t) of local passage, x(0) = x m G E|f, x™ s ^ 



be given. Let a unit tangent vector 5 m to x m G S in with respect to the metric (4-21), and 



its projection 5(0) be given. Then the evolution of S under the linearised flow is estimated 
by 

\S ss (t)\ < exp(-a 2 /isst)C||< s || 2 -\ 



IM*)I < exp(-a 2 ^)C||< n s ||2 1 |a; 



\S c (t) - 5, 



< eC, 



all along the local passage, t G [0, t loc ] 
e < e . 



The constant C is independent of x m , 6 m , t, and 



Proof. Assume 



\5 ss (r)\ < 2C\\xZ 



s II 2 ' 



\S s t(r)\ < 2C, \5 c (t)\ < 2, 



(4.29) 



on r G [0, t\. This assumption holds for small t due the estimates of the initial values in 
lemma 14.51 The constant C is taken from the lemma. 



Then the estimates of lemmata 4.6 and 4.2 yield 



|<Wt)|' < -a/i ss |5 ss (r)|+£Cexp(-a/i ss r)(2 + 2iVC) 
< -afi ss \5 ss (r)\ + eC exp(-afi ss r), 



14 



< 
< 



-afi se \5 se {r)\ + eC exp{-afi se T)\x^\ (2 + 2C||< n s || 2 1 + 2{N - l)C) 
-an se \S se (T)\ +eC , exp(-a / u S £r)||x™ s || 2 " 1 |a;^|. 



The last inequality uses 1 1 ic^ s 1 1 2 < \/N. The new constant is bounded by C < (2 
2NC)\/ r NC and will again be denoted by C. 

Then we can integrate F] the above estimates to obtain 



\8ss(t)\ < exp(-a/i ss t)(|5 ss 
\S se (t)\ < exp(-a/i s ^)(|<^ 



+ eCt) 

+ eCt\\x™\\ 2 



^his can also be seen as a Gronwall estimate for exp(a^, ss t)S ss (t) and exp(a/i s et)S s i(t). 
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For small enough Eq this yields 



|<W*)| < exp(-«Vsst)(|4s(0)| + l) 

< exp(-«Vs S t)(C||xi 1 n s ||2 1 + l) 



< exp(-o; 2 /x ss t)(C + l)\\x 



u,s 1 1 2 ' 



< exp(-aVs^)(C + 1)|K-||2 k 



1 Uin I 



In particular, assumptions (4.29) and the first two claims hold as long as |5 c (r)| < 2, if 
the original constant was chosen larger than 1. 

We use the new estimates of 8 8B , 5 s e, the assumption on 8 C , and the bound on the 



trajectory given by lemma 4.2 to estimate the centre component: 



S c {t)\ < sC l\x ss \ + J2f =1 \x s e\)\x u \\5 c \+eC\x n \[\5 ss \ + J2 1 e=i\^ 



=1 \ x st\ 

+eC{C + l)||< n s || 2 1 \x u \ (exp(-a 2 /i ss r) +J2e=i ex P{-°' 2 f J 'B£ T )\ x ^\ 
< leC (exp(-(a/i ss - /j u )t)\x™\ + Y%=i exp(— af^T") |a^5| 



+£:C(C + l)||x^ s || 2 1 (exp(-(aVss - /i u K)K n | + J2i=i exp(-a 2 /i s ^ 



x 



s£ 



< eC ( exp(-(a 2 /iss - A*u)t) + X)iLi exp(-a 2 /i s ^r) 



with new constant C* < C(C + 3). Thus |5 c (r)| decays exponentially for a close enough 
to 1. Integration yields 

' I4(r)| < 5(7, 



with C < C((a 2 /i ss — /i u ) 1 + X^i^ 2 /^) x )- If is chosen small enough, this shows 



that the assumption (4.29) indeed holds all along the passage and the claimed estimates 
are valid. M 



Theorem 4.8 (local Lipschitz map) The local passage W oc : E 1 " — > S° ut is Lipschitz 



continuous with respect to the metrics (4-21, 4 



There exist (3 > 0, e > and C > such that for all e < Eq the following estimates 
hold for all x m , x m with < x 1 * 1 < x™: 

dist.(z™*,x° u *) < |xjj , pCdist.(x in ,x in ) J 

l^out _ x out) _ (~m _ x in)| < e Q di S t*(x in , X in ). 
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The domain E* n , as defined in (4-25), is given by the local section without the invari- 
ant singular boundaries but including the line (0,x c ) representing the cap of heteroclinic 
orbits. 

The drift in the centre direction can be made arbitrarily small by choosing a suffi- 
ciently small local neighbourhood. The contraction in the transverse directions is arbi- 
trarily strong by restricting the in-section to the part close to the primary object, i.e. the 
stable manifold of the origin. 



Proof. This is a corollary of lemma 4.7 by applying the passage time (4.20) and the 



metric (4.22). Extension to the line (0,x c ) is given by corollary 4.4 



Indeed, in the out-section, the estimates of lemma 4.7, read 





im 


= IM* loc )l 


< 


|j/ u 


o 2 /i ss //i u /^|| ™in II -1 
°ll A u,sll2 ) 




iw 


= IM* loc )l 


< 


|™,in 
l^u 


a 2 [i Be /n u s-i\\ ™in II— l|_.in| 
° Il x u,sll2 l x s^l 


\6r 


-ci 


= i<5 c (t ioc )-ci 


< 


eC. 





With respect to the modified metric (|4.22|) we find using the estimates of lemma 4.2 

r out || 
ss,s II 2 . - out 



\x, 



out I 1 ss 



1*8 



v 2 / N 

i < i + 



in|2(oi/Lt s ^//Lt u — fiss/afiu) |j.™| 2 | IA out l 2 



sC I 



N 



in|2(£+a/W/4 u -/Wa£t u ) I „in|2 



< c- 



N 



\X 



in|2(a 2 /^ ss //i u -5) 



-in 1 2 
x s£l 



< \x™\ 2p C. 

In the second inequality, we introduced a parameter £. The last inequality then needs 

< a 2 /iss//i u -^ =: 0, 

1 < e, 

< ^ + afi se /fi u - fi ss /afx u , E=1,...,N. 
For a close to 1, a suitable £ exists. In fact we can obtain arbitrary 

< < min{ /i ss /> u - 1, Het/pn }. 
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Similarly we find for 1 < £ < N, 



I ™out II \ 2 

I x ss,sll2 | ^out I 



\ X si 



out I \°s£ 



N 



<- I I in|2(a(j ss //i„-fi s( /a(i u ) , \ r | in|2(a/i s| -/^ u -^ s £/a^ u )| in|2 I l^sl 

— I I U I ' / j I U I I s £\ I I i] 



out 1 2 



™m|2 



=1 

N 



| a .in|2(4+aA*ss//iu-/i»{/aift u ) _j_ | a .in|2(^+ttM B f/Mu-rt9l/a/*u)|j.in|2 



1=1 

TV 

ci 2 + £ 



<■ n ^1 Uin|2(a 2 Ms£/Mu-0 

^ U Ar l X u I 



I n I ^ ^S-^l 



< I4 n | 2 ^c. 

This time we need for the last inequality 

< a 2 u. s£ /u. u -£ =: (3, 

1 < £ + a/j ss /fi u - ii st /aix n , 

< £ + a/x^-//x u -Use/an^, £ = 1,...,N. 

Again, for a close to 1, a suitable £ exist. In fact, we can again obtain arbitrary 

< (3 < min{ /x ss //x u - 1, /i^//x u }. 

Now, take a geodesic curve in E m that defines dist*(x m , x in ). The image of this curve 
under the passage \l/ lo€ provides an upper bound on dist„,(x out , x ont ). In both sections the 
x c -component can be separated. Therefore the above estimates on the evolution of the 
tangent vectors immediately yield the claims of the theorem. ixi 



Remark 4.9 In theorem 4-8 , the constant C only depends on the C l bounds on the non- 
linear part of the vector field and the derivatives of the eigenvalues of the linearisation 
along the line of equilibria. The exponent (3 only depends on the spectral gaps. In fact, it 
can be taken arbitrarily in the interval 

< (3 < min{ /x ss //x u - 1, /x s ^//x u }, 

by choosing Eq small enough. 



The last remark provides uniform Lipschitz estimates for the local passages near the 
Kasner circle in Bianchi models, provided they keep a uniform distance from the Taub 
points at which the spectral gap shrinks to zero. 
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In this section we define a global excursion map for trajectories near a primary heteroclinic 
orbit to the Kasner circle, that is from the out section of a local passage to the in section of 
another local passage, both local passages as discussed in the previous section. Combining 
local passage and global excursion we obtain a return map from one in-section to the next, 

* fe := ttf b ottjF : E£ (5.1) 



see figure 5.1 The given heteroclinic orbit corresponds to a fixed origin of this map. 

We prove uniform Lipschitz- and cone properties of the return map, independently of 
the given heteroclinic orbit, as long as the orbit keeps a uniform distance from the Taub 
points. In fact, we prove uniform Lipschitz- and cone properties of the return map under 



the conditions 3.3 on the global excursion. 

This yields a sequence of return maps, with uniform estimates, to every sequence of 
heteroclinic orbits to the Kasner circle that does not accumulate to any Taub point and 



satisfies the local conditions |3.2| at every equilibrium. 

Due to their cone properties, the return maps induce a contracting map on a suitable 
space of sequences of Lipschitz curves. The fixed point provided by the contraction 
mapping theorem then yields the stable manifold of the heteroclinic sequence as claimed 



in theorems 3.1, 3.4 
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Take a sequence p k , k G IN of equilibria on the Kasner circle, not accumulating 



at any Taub point and connected by heteroclinic orbits h (t), lim^oo h (t) = p = 
limt-^-oo h (t) , as in theorem 3.4 Assume that the local conditions 3.2 hold uniformly 
at all p k , in particular sup fcg]N (i u {p k ) / fi ss {p k ) < 1 and mf ke ^ /i se (p k ) > 0. In the Bianchi 



VIo system (2.1) with magnetic field, these conditions are satisfied exactly for a chain of 
heteroclinic orbits not accumulating at Taub points and not containing heteroclinic orbits 



of the magnetic family to equilibria in the intervals /C2, /C5, see figure 2T In particular, 
the conditions hold for the period 3 cycle. 

The previous section then applies to all p k and the coefficients £q,(3,C of the local 



estimates of theorem 4.8 can be taken independent of k, see remark 4.9 



Note the order of fixing the rescaling parameters: First Eo resp. e is fixed small 
enough to yield our estimates of the local passages \I/j? c with small Lipschitz constants, in 
particular eC <C 1 in theorem 4.8 This amounts to a choice of the sections £ in (pfc) and 
E out (pfc) in the original (unsealed) coordinates and also fixes the global excursion maps 
\£f lob 

k 



Due to the non-Euclidean metric used in theorem 4.8, we have to restrict our local 

out 



passage map to ^/ lo< 



^k 
S out 



k 

\pk) 



Ej? -> £° ut by 



\ (•£'u,S; <^c. 
{ (•Z\SS,S ; 



e ^ n {Pk) I 

, G S out (p fe ) 



or x u x sl ■ ■ ■ x sN ^ }, 



(5.2) 



X Q 



= or X ss X s i • • - x sN 7^ }> 
see (4.25). Then a sufficiently small upper bound for xJJ 1 is chosen, i.e. ^/j? are restricted 
to smaller sections 



= { x G S^ 1 I < x u < 5, < x s£ <5, \x c \ < 5 } = ((0, <5] JV+i U 0) x [-5, 5}. (5.3) 

This makes the contraction of the local passage as strong as we like without changing 
. It also ensures that trajectories of interest stay close to the Kasner caps of 



lob 



heteroclinic orbits and therefore the global excursions ^f oh on the domain of interest are 
as close to the Kasner map as we like. It also ensures that all non-singular trajectories in 
these domains indeed return to the following in-sections 1 . 



The global conditions 3.3 hold accordingly: the invariant subspaces, (glob-i), are 
those of the Bianchi system; the uniform bound, (glob-ii), is fixed by choice of a uniform 
size Eq of all local neighbourhoods; and expansion, (glob-iii), is given by the Kasner map. 
Uniform expansion again needs a uniform distance from the Taub points. 
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The following lemma relates the global excursions to the new metric used for the 
local estimates. 



Lemma 5.1 Let the global conditions 3.3 be satisfied for the sequence ^f° h '■ E°, ut — > S in 



fe+i 



of global excursions. Then condition (glob — ii) also holds with respect to the new metrics 



(4.21, 4.22) 



Proof. Due to the invariant subspaces, (glob-iii), the linearisation D^f oh are diagonal 



at the origin and close to diagonal in the neighbourhoods of interest. Moreover, the 

f gi< 

k 



transverse components of ^/f oh have the form 



ss, si, ... , sN, 



with smooth ^f lob . Now note the definitions (4.21, 4.22) of the new metric. The bounds 



(glob-ii), on first and second derivatives of ^f° b , yield uniform bounds on \l/f lob and 
their first derivatives. Thus, the ratio of the coefficients of the metric at an arbitrary 



x e E? to the coefficients at \I/ 



glob . 



k 

ob\— l 



X 



E E 1 ^ is between 



glob\ 



- 1 II — 1 1 



glob I 



and 



|\l/| lob ||, that is between M~ 2 and M 2 , for a uniform constant M. This im- 
mediately yields new uniform bounds on the derivatives of \l> glob with respect to the new 
metric. M 



Now we can proceed along the lines of [LHWGlT] (section 4) to establish the existence 
of stable manifold by a graph-transform approach. 



Lemma 5.2 Assume conditions 5\£ on the local passages and conditions 5\3 on the global 
excursions. 



Then the return maps (5.1) are Lipschitz continuous with respect to the metric (4-21). 
Furthermore, there exist e > 0, 8 > 0, < a < 1, K u s > 1, and K c > (1 — ct 2 ) -1 > 1, 
such that the following cone conditions hold for 



* fc = *f b o tff : E- -> Ef +1 



Here E™ are the in-sections (5.2) corresponding to the choice of e, and E™ are suitable 



subsets of the form (5.3) 
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^u,s Expansion 
of x c 
Contraction 
of x u s 




^ty(x) 



Figure 5.2: Cone properties of the return map 



The cones are defined for x G S m ( omitting the index k) as 



The cone conditions are 



{s G S in | dist* (x us , x us ) < o~\x c x c \\, 
{5; G S in | \x c — x c \ < crdist*(x UiS — x u>s )}. 



(5.4) 



(i) Invariance: *(C£) n S m C (int C&J U {^a;} and * _1 (CJJ^) n £ m C (int Q s ) U {x}; 

(wj Contraction & Expansion: For all x G C£ we /iai>e expansion in the centre direction: 
\(^x) c — (^x) c \ > -K" c |5c — a? c | an d f or <dl ^x G C^ s x we have contraction in the 
transverse directions: dist*(x u , s , x MiS ) > X U)S dist*((^ r x) U)S , (^)«,s)- 



5.5 



T/jq/ hold for all, x, x, tyx, ^x G E m . See also figure 

The coefficients a, 5 only depend on Eq and the uniform expansion (glob-iii), that is 
the distance to the Taub point in the Bianchi system. 



Proof. Lipschitz continuity of the return map follows directly from Lipschitz conti- 
nuity of the local passage ^j? , see theorem 4.8, as the global excursion \I/f lob is smooth. 
To simplify notation, we drop the index k from now on. All estimates will be uniform in 
k. 

The cone conditions require the expansion in x c -direction given by (glob-iii), cor- 
responding to the expansion Kasner circle induced by the Kasner map of the Bianchi 
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system. In fact, (glob-iii) states that we have 



^glob/^out 



0, x s 



out 



0,0 = K n = 0,x- = 0,$(O) 



(5.5) 



with Lipschitz continuous $, $ 1 . The Lipschitz constant of $ 1 is less than L < 1 
independent of k. Note again the invariant boundaries, (loc-v), (glob-i). Therefore we 



can write, as in the proof of lemma 5.1 



# glob (x ss , x s , x c ) = (0, 0, $«)) + * glob (x ss , x 8 , x c )x SSyS 



with a smooth matrix \l/ glob (x ss , x 8 , x c ) and vector x SSjS = (x ss ,x s ). 

Consider now two points x,x G S out . Choose geodesic paths 71 from to x sss and 



72 from x SSjS to x SSjS , both with respect to the new metric (4.22). Then 



dist^ glob (x),^ glob (a;)) 



f-dist„(0,x SSiS ) 1 

+ / — dist,(v]/ glob (7 1 ( s ),5 c ),vl/si°b (7i(s));Ec))ds 

pdist* [X ss S ,X SS s ) J 

/ — dist,(^ glob (7 2 (s),5 c ),^ glob (x))d S 



ds 



< \$>(x c ) - $(x c )\ 

fdist»(0,a;ss,s) ^ 



+ 







ds L dt 



dist* (* glob ( 7l (s) , t) , ^ glob (71 (s) , x c ) ) dt ds 



fdist*(x aa , s ,x SSiS ) j 

- / Q — dist,(vl/ glob ( 72 (s),x c ),vI/ glob (a;))ds, 



and we obtain the following Lipschitz estimate 

I dist*(^ glob (x),^ glob (x)) - |$(z c ) - $(x c )| I 

< C* g (dist* (0, x ss s ) \x c x c \ + dist*(x SSjS , x SSiS )) , 
^ C g (ll^ssjsll |-^c -^c| "I - dist* (x SSiS , x SS)S ) ) , 

with C gloh only depending on the uniform bounds on ||D\l/ glob ||, ||D 2 \I/ glob || with respect to 



the new metric provided by lemma 5.1 The last inequality used the trivial upper bound 



(4.24) on the distance from the origin in the new metric. 

Using \I/ loc (x) and ty loc (x) instead of x and x we get a similar estimate for the return 
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map V(x) = ^ glob (^ loc (x)): 

| dist,(tf (x), - |$(^ loc (5) c ) - $(* loc (x) c )| | 

= | dist*(^ glob (^ loc (x)) - ^ glob (^ loc (x))) - |$(^ loc (x) c ) - $(^ loc (x) c )| | 

< C glob (||^ loc (x) ss , s |||^ loc (x) c - ^ loc (x) c | +dist,(^ loc (x) ss , s ,^ loc (x) ss , s )) (5.6) 

< C glob (\x n \P\\x\\(l + eC)disU(x,x) + \x u \^CdisU(x 7 x)) 

< C Ietuin \x u \Pdist*(x,x). 



The second last inequality uses the estimates of the local passage of corollary 4.4 and 



theorem 4.8 for the choice (w.l.o.g.) < x u < x u . Note that the estimates of theorem 4.8 



are used in the form 

|^ loc (x) c - ^ loc (x) c | < eCdisU(x,x) + \x c -x c \ < (eC + l)dist*(x, x), 
dist*(^ loc (£) SSjS ,^ loc (;r) ss , s ) < \x^Cdist m (x,x). 

The constant (7 rcturn is uniform in x, x in the in-section, and the omitted number k of the 
section along the heteroclinic chain. Because < (3 < min{/i ss //i u — l,/i s ^//i u }, we have 
an arbitrarily strong contraction for x u < 5, if we choose 5 small enough. 

The map $ given by ( 5.5[ ), i.e. the Kasner map in the original Bianchi system, is 



expanding, see condition (glob-iii): 

|$(a)-$(6)| > L _1 |a-6|. 

for some uniform constant L < 1. 

Now choose K c with 1 < K c < L~ l , and o with < a < 1 such that K c (l — a 2 ) > 1. 
(The last relation is needed to obtain a contraction in theorem 5.3[ ) 



Consider the cone in centre direction with opening $ > 0, i.e. dist*(x UiS , x U)S ) < 



i9|:r c — x c \. Then (5.6) using the local Lipschitz estimate of theorem 4.8 yields 



#(x)) c | > |$(^ loc (x) c ) - $(^ loc (s) c )| - C Tetmn \x u \PdisU(x,x) 

> L- 1 |^ loc (x) c - ^ loc (x) c | - C retw:n \x u fdisU{x,x) 

> L~ 1 \x c — x c \ - L-^Cdist^x.x) - C^^lx^^st^x) 

> (L- 1 - (L^eC + C ietwcn \x u \P) (1 + $)) \x c - x c \. 

For £ and 5 chosen small enough, using \x u \ < 5, we can achieve 

K c < L- 1 - (L^eC + C vctuin \x u \P) (1 + l/a), 



(5.7) 
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yielding the expansion not only in the cone with $ = cr < 1, but also outside the cone 
Q' s , with i? = 1/cr. 



Furthermore, using again (5.6), we see the invariance of the cones. Indeed, assume 



again dist*(x UiS , x U)S ) < $\x c — x c \, then we have 

dist,(tt(z) U)8> tt(a;) U)B ) < C^ m \x u fdisU(x,x) 

< C return |x u |^(l + i?)|x c -x c | 

< C return |x u | /3 (1 + fyK-^ix),. - *(x) c |. 

The last inequality uses the expansion in x c , thus it is valid for *& < 1/cr. We choose 8 
small enough such that C xetwn \x u \" K~ l < cr/(l + cr). Due to the monotone increase of 
#/(l + •&) we also have C IetuTn \x u f K' 1 < 0/(1 + ti) for all & > a. Thus we obtain the 
cone invariance 

dist*(*(x) UiS) #(x) U)S ) < 0|tt(z) c - (z) c | (5.8) 

for all (J < I? < 1/(7. 

The choice $ = a yields (forward) invariance of the cone C% and the choice d = 1 jo 
yields (backward) invariance of the cone C^. Note that the cone invariances are in fact 
strict as claimed in the lemma. The above estimates are strict inequalities for x ^ x. 

Now consider the cone in transverse direction, that is ty(x) G C^, which amounts to 
|\E'(x) c — x I / (a;) c | < crdist*(^ r (5) U)S , \&(x) U)S ). We have already established invariance. Thus 



|x c — x c \ < cr||s UiS — x us \\ and estimate (5.6) yields 



dist*(#(x) U)S ,#(x) U)S ) < C rctum |x u |^dist,(x,a;) 

< C rcturn |x M |^(l + Cr)dist*(i u , s , X u , s ). 

This is the claimed contraction, K~l = C ietuIIl 8"(l + cr) < 1, for 8 small enough. ix 



Theorem 5.3 Assume conditions 3.2 on the local passages and conditions 3.3 on the 
global excursions. 

The (local) stable set of the origin under the sequence of return maps is given by 

~\ Alloc — / / k k k\ \ k — k / k k \ \ 
yv k — U x u? x s£i x c) I x c — x cl x u? x s£JJ ■ 



The functions x^ are Lipschitz continuous with respect to the metric (4-21). Furthermore, 
z*(0) = and #(W£ oc ) C W^. 
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Proof. The idea of the proof is to define a graph transformation on the space of sequences 
of Lipschitz-continuous graphs {x UjS H- x c = Cfc(^u,s)) I k G IN} by the inverse return maps 
^l/^ 1 . The uniform cone invariance provided by the previous lemma will ensure that the 
Lipschitz property of the graphs is preserved. Due to the expansion/contraction conditions 
of the previous lemma, the graph transformation turns out to be a contraction on the 
space of sequences of Lipschitz-continuous graphs. The fixed point of this contraction 
then yields the claim. 

To make this idea precise, consider the Banach space of Lipschitz-continuous func- 
tions 

X = { ( : (0,5] 1+N U0^ [-6,6], x u , s = {x u ,x s ) ^x c = C(x u ,s) 

such that Lip(C) < o and £(0) = } 



with sup-norm. The parameters 8, a < 1 correspond to those of lemma 5.2 Lipschitz 



continuity is considered with respect to the metric dist* given by (4.21). Consider also 
the space of sequences 

X M = { (Cfc)fceM I Ck G X } 

with sup-norm. 

Define maps Gk ■ X — > X as graph(GfcCfc+i) : = ^ 1 g r aph(Cfc+i), i-e. as the transfor- 
mations of the graphs of the functions in X. More precisely 

GfeCffeH^.CWJluJ : = (^Ks.CW)),, forx UiS ^0, 

GC(0) := 0. 

The first equation implicitly assumes that (x UjS , C(^u,s)) has a pre-image under \1> and that 
it lies in the domain. The second equation just gives the pre-image of the origin under 
Note the restriction to non-negative x u , x 8 i, o; s n consistent with the invariant boundaries 
(loc-v), (glob-i). 

We will prove the following claims, uniformly in the index k, (which is dropped from 
now on to simplify notation) 



(i) domain of definition: for all ( G X and x u>s G (0, 5] N+1 there exists x u>s G (0, 5} N+1 , 
such that (* -1 (£u, s ,C(^u,s))) U!S = %u,s- 

(ii) well-definedness: for all ( G X and x UjS ,x UiS G (0,5} N+1 the following holds. If 

(^ _1 (^u,s,C(^u,s))) us = (^ _1 (^u,s,C(5u,s))) us G (O,^^ 1 then already x U:S = x u>s . 
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Conditions (i) and (ii) yield a well denned function G( with (G()(0) = for every ( G X. 

(iii) Lipschitz property: for all £ G X the function G( is again Lipschitz continuous 
with Lipschitz constant Lip(G() < o. Note that the Lipschitz property is again 
considered with respect to the metric dist*. 

(iv) contraction: The exists a constant < K < 1 such that for all (, ( G X the estimate 
IIGC-GCIIsup^llC-Cllsup holds - 



Conditions (i)-(iii) prove that the graph transformation G indeed maps Lipschitz contin- 
uous functions in X to Lipschitz continuous functions in X, with respect to the metric 
dist*. Condition (iv) provides a contraction. Uniformity of bounds yield a contraction 
on the space X w of sequences. If all four conditions hold, then by contraction-mapping 
theorem there is a unique fixed point, i.e. a sequence of Lipschitz continuous function 
C fc * with G k Q +1 = Q. 

Its graphs form a forward invariant set under \l/ composed of local manifold. It is 



also the stable set of the origin due to the cone conditions of lemma This yields the 
claim of the theorem. Therefore it remains to prove (i)-(iv): 

(i) Let (6l and x UjS € (0, 5} N+1 be given. The straight line {x UjS } x [—5, S] is 
contained in the cone C , ( C Xus0 )- We use lemma 5.2 \l/(x UiS ,0) G S m by invariance and 
contraction of the cone Cq' s . Thus, by invariance and Expansion of C? x ^, the image 
of the straight line {x u>s } x [—5,5] under \1/ contains a curve in C^,^ Xus0 ^ connecting the 
extremal planes {x c = ±5}. By the intermediate value theorem this curve must intersect 
the graph of (. 

(ii) Let C € X and x UjS ,x UjS G (0,5] N+1 be given with (* _1 (x U)S , C(^u,s))) UiS = 

(*- i (£u,s,c(5u, s ))) U)8 G(o,5r+ i . 

Then ^ _1 (£ u ,s, C(5 u ,s)) G Q,-i (Xu siC(a;u s)) , and by cone invariance (z UiS , C(^u,s)) G 
C? s . The Lipschitz-bound on £ G X on the other hand implies (x U}S , C(^u,s)) £ 

^(Ju, s ,((i u , s ))) thus (x u ,s, C(^u,s)) — (^u,sj C( x u,s))- 

(iii) Again, the Lipschitz-bound on ( G X translates to (£ u ,s, C(^u,s)) £ ^(X s C(^u s )) 



for all Cone invariance and lemma 5^2, immediately yield the Lipschitz bound on 

G(. 
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(iv) The origin is fixed by construction, thus we only have to estimate the dis- 
tance of the nonsingular part. Let (,( G 1 and x UjS ,x U)S G (0,5] N+1 be given with 

(^- 1 (x u , s ,C(x u , s ))) u , s = (^- 1 (x u , s ,C(xu,s)))u, s e (0,5] N+1 . 

Again, this implies l I /_1 (a; UiS , C(x u ,s)) G C*-i(i us ((i us ))i anc ^ by cone invariance we 
have (x UiS , C(^u,s)) C(L US C(^ us ))' Thus we can estimate 

|C(f u ,s) -CWI < IIC - Cllsup + o-dist*( 

■^U,S) -£u,s ) 

< IIC- Cllsup + a 2 |C(£ U)S )-C(^,s)| 

The first inequality uses the Lipschitz bound on ( G X whereas the second one uses the 
aforementioned cone C? r , y,. We obtain 



|C>u,s)-CKs)| < IIC- Cllsup- 

On the other hand, the expansion of C , ^_ 1 , a . us f^^) under \l> yields 



(gc-gc) ((r 1 (. ¥1 ((i u , 5 ))) u 

(^- 1 (x u , s ,C(x u , s ))) c - (^- 1 (x UiS ,C(5 u ,s)) 

C("^u,s) C("^u,s) 



< 



< 



K, 



^iic-cii 



sup- 



Lemma 5.2 provides constants K c , a with iT c (l — a 2 ) > 1. Therefore the last estimates 



yield the claimed contraction, k = 1/(K C (1 — a 2 )), and this finishes the proof. 



tX] 



With theorem 5.3 we have finally proved the main theorems 3.1, 3.4 as formulated 
in section [3] 



6 Discussion and outlook 

Unfortunately the metric used to obtain the contraction in the proof of the main theorem 
of this paper is singular on the invariant subspaces. We have no result on the way in 
which the manifolds constructed approach these boundaries. The only exception is the 
heteroclinic cycle itself. The cap of heteroclinic orbits corresponds to the line {0, x c } and 
the new metric is regular there. In fact, the manifold is nicely attached to the given 



heteroclinic chain. Moreover, in the proof of theorem 5.3 we could restrict to very small 
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neighbourhoods of the primary heteroclinic chain, i.e. 5 — > 0. Then we can choose 
arbitrarily small Lipschitz bounds on the functions considered, i.e. o — > 0. Thus the 
manifolds constructed are tangent to the fibre {x c = constant} at the heteroclinic chain. 

For completeness the following subtlety should be mentioned. The set of points in the 
domain of definition of the dynamical system corresponding to Bianchi type IX vacuum 
solutions or Bianchi type VIo solutions with magnetic field is an open subset bounded 
by invariant manifolds and it lies on only one side of these manifolds. The fact that the 
global excursion map has its image on the correct side of these manifolds is not mentioned 
in the analytical treatment above. Nevertheless it follows immediately from the nature of 
the underlying geometrical problem. 

Up to this point vacuum models of type IX were replaced by Einstein-Maxwell models 
of type VIo and one non- vanishing magnetic field component. Now some generalizations 
will be mentioned. In 1.1 1\\ (i 1 lj some results were obtained for type IX solutions with 
perfect fluids having a linear equation of state p = (7 — l)p. Restrictions had to be 
imposed on the value of 7. The techniques developed in this paper allow these results 
to be generalized to cases where these restrictions are relaxed. The four-dimensional 
dynamical system is replaced by a five-dimensional one and at each Kasner point there 
is an additional eigenvalue 3(2 — 7). This situation can be treated for all 7 < 2. In 
particular the method applies for all values of 7 in the physical range [1, 2] except for the 
case 7 = 2 where the dynamics is known to be very different. Bianchi type VIo solutions 
with a perfect fluid and a magnetic field can be treated in a very similar way. The 
additional eigenvalue arising from the fluid is the same as in the case without magnetic 
field jLKW95j . 

It is possible to formulate the Bianchi type II models with a magnetic field as a 
five-dimensional dynamical system [LeB97j. In this approach the magnetic field has only 
one non-zero component in the frame used but the metric has a non-zero off-diagonal 
component in that frame. The eigenvalues of the linearisation about a Kasner solution 
are given by 

3pi, Qp 2 , 3(p 3 -p x ). (6.1) 

An important qualitative difference to the model of type VIo is that for some regions of 
the Kasner circle the stable manifold of the Kasner solution is two-dimensional. Thus 
in general the results of this paper do not apply to heteroclinic chains for the Bianchi 
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type II model with magnetic field. In a similar way it is possible to formulate the type 
I models with a magnetic field as a five-dimensional dynamical system with only one 
component of the magnetic field being non-zero |LeB97] . In this case all the off-diagonal 
metric components are non-zero in general. Again it happens that the stable manifold 
can be two-dimensional. Note that it has been shown in |LKW95] that it is not possible 
to have solutions of the Einstein-Maxwell equations of Bianchi type VIII or IX with a 
non- vanishing pure magnetic field. 

Up to now there is no generalization of the results of |LHWGlT] to oscillatory models 
of Bianchi class B. In fact it would be very interesting to have such results for Bianchi 
type VT_i where oscillatory solutions are expected to exist. One obstacle is the existence 
of stable manifolds of dimension greater than one as in the examples with magnetic field 
above. Another is that invariant manifolds of the type which played such an important 
role in the proofs of this paper do not appear to exist for models of Bianchi class B. 

In the case of Bianchi type IX vacuum models it has been proved that the a-limit set 
of each solution belongs to the union of points of type I and type II |Rin01j . Interestingly it 
is not known if the corresponding statement holds for the superficially similar type VIII. 
This contrasts with the fact that the results for type IX in (LHWGllj extend almost 
without change to type VIII. It is easy to formulate an analogue of the result of |Rin01] 
for solutions of type VIo with magnetic field and it would be interesting to investigate 
whether it holds, especially since this might throw some new light on the unsolved Bianchi 
VIII problem. 

To sum up, it is clear that the above complex of problems represents a promising 
opportunity to learn about the related questions of the BKL conjecture, the dynamics of 
Bianchi models near the initial singularity and the stability of heteroclinic cycles in more 
general dynamical systems. 
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